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Abstract 

Background: The eyes and skin are obvious retinoid target organs. Vitamin A deficiency causes night blindness and 
retinoids are widely used to treat acne and psoriasis. However, more than 90% of total body retinol is stored in liver 
stellate cells. In addition, hepatocytes produce the largest amount of retinol binding protein and cellular retinoic 
acid binding protein to mobilize retinol from the hepatic storage pool and deliver retinol to its receptors, 
respectively. Furthermore, hepatocytes express the highest amount of retinoid x receptor alpha (RXRa) among all 
the cell types. Surprisingly, the function of endogenous retinoids in the liver has received very little attention. 

Results: Based on the data generated from chromatin immunoprecipitation followed by sequencing, the global 
DNA binding of transcription factors including retinoid x receptor a (RXRa) along with its partners i.e. retinoic acid 
receptor a (RARa), pregnane x receptor (PXR), liver x receptor (LXR), farnesoid x receptor (FXR), and peroxisome 
proliferator-activated receptor a (PPARa) has been established. Based on the binding, functional annotation 
illustrated the role of those receptors in regulating hepatic lipid homeostasis. To correlate the DNA binding data 
with gene expression data, the expression patterns of 576 genes that regulate lipid homeostasis were studied in 
wild type and liver RXRct-null mice treated with and without RA. The data showed that RA treatment and 
RXRa-deficiency had opposite effects in regulating lipid homeostasis. A subset of genes (114), which could clearly 
differentiate the effect of ligand treatment and receptor deficiency, were selected for further functional analysis. The 
expression data suggested that RA treatment could produce unsaturated fatty acids and induce triglyceride 
breakdown, bile acid secretion, lipolysis, and retinoids elimination. In contrast, RXRa deficiency might induce the 
synthesis of saturated fatty acids, triglyceride, cholesterol, bile acids, and retinoids. In addition, DNA binding data 
indicated extensive cross-talk among RARa, PXR, LXR, FXR, and PPARa in regulating those RA/RXRa-dependent 
gene expression levels. Moreover, RA reduced serum cholesterol, triglyceride, and bile acid levels in mice. 

Conclusions: We have characterized the role of hepatic RA for the first time. Hepatic RA mediated through RXRa 
and its partners regulates lipid homeostasis. 
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Background 

The eyes and skin are obvious retinoid target organs. Vita- 
min A deficiency causes night blindness and retinoids are 
widely used to treat acne and psoriasis. However, more than 
90% of total body retinol (retinylpalmitate, the storage form) 
is stored in liver stellate cells [1]. In addition, hepatocytes 
produce the largest amount of retinol binding protein and 
cellular retinoic acid binding protein to mobilize retinol 
from the hepatic storage pool and deliver retinol to its re- 
ceptors, respectively [2]. Furthermore, hepatocytes express 
the highest amount of retinoid x receptor alpha (RXRa) 
among all the cell types. Surprisingly, the function of en- 
dogenous retinoids in the liver has received very little atten- 
tion. Thus, the current study aims to identify the bona fide 
RXRa and RARa targets in the liver. 

The broad and complicated roles of retinoids can be 
explained by the presence of multiple receptors for retin- 
oic acid (RA), the biological active form of retinol. The re- 
ceptors for RA are retinoic acid receptor (RAR) as well as 
retinoid x receptor (RXR) [2]. In addition, RXR is essential 
for many other receptors to function. These receptors for 
RA belong to a nuclear receptor family whose members 
are transcriptional factors. Thus, RA exerts its biological 
effects by regulating gene expression. RXR is unique in 
that it not only forms homodimers, but also dimerizes 
with other nuclear receptors, which include receptors for 
fatty acids (peroxisomal proliferator-activated receptors, 
PPARs), bile acids (farnesoid x receptor, FXR), oxysterols 
(liver x receptor, LXR), xenobiotics (pregnane x receptor, 
PXR, and constitutive androstane receptor, CAR), vitamin 
D (vitamin D receptor, VDR), and RA (RAR). Hence, most 
RXR partners participate in regulating lipid homeostasis. 
Within these heterodimers, RXR can be either a permis- 
sive or a silent partner. When RXR serves as a silent part- 
ner, the heterodimer does not respond to RA. When it is a 
permissive (active) partner, RA and the ligand for the 
heterodimeric partner can both activate the heterodimer. 
For example, RXR is a permissive partner for PPAR [3]. 
Similarly, heterodimeric complexes of RXR with LXR [4] 
or FXR [5] also retain RA responsiveness. Furthermore, 
retinoids also activate PXR, VDR, and CAR thus are able 
to thus regulate xenobiotic metabolism and potentially 
their own oxidation [6-8]. Since most of these receptors 
are abundandy expressed in the liver, the endogenous RA 
may regulate many hepatic nuclear receptor-mediated 
pathways. Therefore, the role of RA in the liver is unpre- 
dictable. In order to understand the endogenous function 
of RA and its receptors, it is crucial to identify RA recep- 
tor targets (genes and pathways) genome-wide. 

RXRa is highly expressed in the liver [9]. Liver specific 
RXRa-deficient mice have increased serum triglyceride and 
cholesterol levels [10,11]. In addition, lack of hepatic RXRa 
increases sensitivity to alcohol- and non-alcohol-induced 
steatosis and steatohepatitis [12,13]. Besides regulating 



lipid metabolism, hepatocyte RXRa also controls xeno- 
biotic [14-16], carbohydrate [17], and amino acid metabol- 
ism [17]. These findings indicate that RXRa-mediated 
signaling has a huge impact on maintaining liver health 
and in regulating many disease processes. 

To understand the global roles of RXRa and RARa at 
the genomic level, chromatin immunoprecipitation using 
anti-RXRa and -RARa antibodies followed by sequencing 
(ChlP-seq) was performed. Since RXRa is an essential part- 
ner for other nuclear receptors, we compared ChlP-seq 
data to RXRa binding locations with locations from previ- 
ous studies for PXR [18], LXR [19], FXR [20], and PPARa 
[19]. Meanwhile, the expression levels of the genes respon- 
sible for lipid homeostasis were studied in wild type and 
hepatic RXRa-deficient mouse livers. Both genome-wide 
DNA-binding and hepatic gene expression data were used 
to define the role of RA in the liver. Our data uncovered 
the unknown function of retinoic acid and RXR vs. RAR in 
the liver. Using different approaches, we showed for the 
first time that retinoic acid-activated RXRa and RARa have 
distinct effects. Moreover, the action of retinoic acid in the 
liver is to regulate lipid homeostasis specifically by redu- 
cing serum cholesterol, triglyceride and bile acid levels. 
The data provided may lead to future development of syn- 
thetic retinoid that can target metabolic syndrome or other 
types of lipid-associated health issues. 

Results 

Genome-wide binding of RXRa, RARa, PXR, LXR, FXR, and 
PPARa in mouse livers 

To understand the global roles of RXRa and RARa at 
the hepatic genome level, ChlP-seq was performed using 
anti-RXRa and -RARa antibodies. Single read sequen- 
cing yielded 18 and 32 million uniquely mapped reads 
for RXRa and RARa, respectively. After filtering by in- 
cluding peak scores that were greater than 20 and dis- 
tance within 10 kb from the transcriptional start site, 
17,973 peaks were detected for RXRa and 18,697 peaks 
for RARa. Since RXRa is an essential partner for other 
nuclear receptors, we compared ChlP-seq data of RXRa 
with those of PXR [18], LXR [19], FXR [20], and PPARa 
[19]. Our data showed that the numbers of peaks, which 
were commonly bound by RXRa and other nuclear re- 
ceptors, were 6,577 for RARa, followed by 5,154 for 
PPARa, 2,846 for FXR, 1,190 for LXR, and 868 for PXR 
in the mouse liver genome. RARa had less than 50% 
overlapping bindings with RXRa while PPARa, FXR, 
LXR, and PXR had over 85% overlapping bindings with 
RXRa (Figure 1A). These findings suggest that RXRa is 
indispensable for the function of PPARa, FXR, LXR and 
PXR in the mouse liver. In contrast, RARa may work in 
the absence of RXRa. More than 4000 of RARa binding 
genes were bound by RXRa, and occupied 56% of total 
RXRa bindings, followed by PPARa (43%), FXR (25%), 
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Figure 1 Histograms showing genome-wide binding of RXRa/nuclear receptor in wild type mouse livers. (A) The overlapping binding 
between indicated nuclear receptor (NR) and RXRa. (B) The preferential dimerization of RXRa with nuclear receptors. The number and percentage 
of genes bound by RXRa and indicated nuclear receptors are shown in black. The binding occurs in the same location within a gene. I: the 
number and percentage of overlapping binding between RXRa and any one of the indicated nuclear receptors. II: the number and percentage of 
overlapping binding between RXRa and any one of PXR, LXR, FXR, and PPARa. 



LXR (12%), and PXR (8%). The well-known lipid regulators 
including PPARa, FXR, LXR, and PXR bound to almost 
50% of RXRoc-bound regions suggesting the extensive role 
of RXRa in regulating lipids. 

Cluster and principal component analysis (PCA) were 
performed using the binding data. Overlapping bindings be- 
tween RXRa and RARa, PPARa, FXR, LXR, and PXR were 
analyzed. According to the length of the vertical branches 
(Figure 2A), LXR and PXR were clustered into one basic 



group. FXR was clustered into a subgroup with LXR 
and PXR. PPARa, PXR, LXR, and FXR were classified 
into another group. However, RARa was distant from 
the others. The same datasets were subjected to PCA 
analysis. After dimension deduction, two components 
were picked to describe the global properties of the 
samples. The two-dimension score plot showed that 
LXR and PXR had the shortest distance, implying their 
properties on the genome-wide binding level were most 
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Figure 2 Genome-wide profiling of RXRa/nuclear receptor binding sites. (A) A dendogram was generated by cluster analysis to show RXRa/ 
nuclear receptor (NR) binding sites. Each horizontal line represents a binding site. White regions indicate there is no peak called. The lengths of 
the vertical lines represent the distance between the sets of target genes for the nuclear receptors when clustered. (B) A score plot of PCA 
analysis for RXRa/nuclear receptors binding sites. Scores of component 1 and 2 were obtained from the linear combination of all the nuclear 
receptors binding sites. 
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similar (Figure 2B). The distance between FXR and LXR/ 
PXR was shorter than the distances between other recep- 
tors and LXR/PXR. RARa was far from all other nuclear 
receptors. Thus, the information generated from the PCA 
and cluster analysis was consistent. The difference among 
the binding profiles of PXR, LXR, FXR, and PPARa were 
described by component 2 while the difference between 
the binding profile of RARa and those four nuclear recep- 
tors were described by component 1 in the PCA. 

Function analysis of the genes that have overlapping 
binding sites between RXRct and each of RARa, PXR, LXR, 
FXR, and PPARa 

To understand the potential biological role of genes likely 
to be targeted by RXRa-heterodimers, DAVID functional 
annotation was performed. The top ten processes for each 
of the heterodimers formed by RXRa and other nuclear 



receptors were investigated in this study (Table 1). Al- 
though common pathways among these heterodimers, such 
as carboxylic acid catabolic process, organic acid catabolic 
process, and oxidation reduction, appeared in the top 10 
pathways, the heterodimers still have their unique biological 
functions. RXRoc-RARoc uniquely regulated protein trans- 
port, localization, and RNA processing. The heterodimers 
of RXRa and PXR, LXR, FXR, and PPARa have a great im- 
pact on various lipid processing pathways. For example, 
RXRa-PPARa and RXRa-LXR tend to bind to the genes in- 
volved in acylglycerol, glycerol ether, and neutral lipid meta- 
bolic processes, which are fatty acid derived. In addition, 
RXRa-PPARa distinctly binds to the genes that are in- 
volved in fatty acid metabolic processes, and RXRa-LXR is 
more prone to bind to the genes that participate in steroid 
metabolic process. In addition to steroid metabolic process, 
RXRa-FXR also tends to bind to genes that dictate lipid 



Table 1 Common and unique biological function of genes bound by RXRa coupled with RARa, PXR, LXR, FXR and PPARa 

Pathways RARa PXR LXR FXR PPARa 

Acylglycerol metabolic process 1.9E-07 3.8E-10 

Neutral lipid metabolic process 2.7E-07 94E-10 

Glycerol ether metabolic process 2.7E-07 1E-10 

Organic ether metabolic process 4.5E-1 1 

Carboxylic acid catabolic process 6.1E-11 7.3E-07 2.8E-08 1.9E-09 6.1E-11 

Organic acid catabolic process 6.1E-11 7.3E-07 2.8E-08 1.9E-09 6.1E-11 

Coenzyme metabolic process 9.4E-13 2E-07 1.7E-10 

Cofactor metabolic process 1.6E-17 3.0E-06 2.1E-07 2E-1 1 

Fatty acid metabolic process 1.8E-12 

Oxidation reduction 1.7E-23 1E-11 1.4E-16 6.1E-14 34E-28 

Sterol homeostasis 7.2E-08 

Cholesterol homeostasis 7.2E-08 

Lipid homeostasis 4.8E-07 

Lipid transport 2E-07 

Hexose metabolic process 2.1E-07 

Steroid metabolic process 2.6E-1 1 2E-08 

Monosaccharide metabolic process 14E-07 

Pyruvate metabolic process 1.3E-05 

Acute-phase response 3.2E-05 

Cellular carbohydrate biosynthetic process 5.9E-06 

Acute inflammatory response 7.3E-07 

Glucose metabolic process 0.00001 

Response to wounding 3.5E-07 

Electron transport chain 2.7E-10 

ncRNA metabolic process 4.2E-15 1.8E-09 

mRNA processing 6.7E-16 

Translation 7E-16 

RNA processing 8.8E-16 
Numbers in the table represent the p value given by DAVID. 
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transport and regulate steroid and monosaccharide homeo- 
stasis. RXRoc-PXR targets specific biological processes like 
regulating pyruvate and carbohydrate biotransformation, 
but it also has a role in acute-phase response, acute inflam- 
matory response, and response to wounding. Taken to- 
gether, RXRa-RARa targeted protein and RNA processes 
while RXRa-PXR/LXR/FXR/PPARa targeted homeostasis 
regulation of small molecules, which include monosacchar- 
ide and lipids. 

Global profiling of the expression of lipid homeostasis 
genes in wild type and RXRct KO mice treated with and 
without RA 

Hepatic RXRa KO mice have elevated serum cholesterol 
and triglyceride levels [21], and the majority of PXR, LXR, 
FXR, PPARa-bound genes were also bound by RXRa. We 
next tested a hypothesis that RA and hepatic RXRa could 
regulate lipid homeostasis in the mouse liver via RXRa 
and its heterodimeric partners. The expression levels of 
the lipid homeostasis genes (579) in the KEGG pathway 
database were studied in wild type and hepatic RXRa KO 
mice treated with and without RA. The PCA score plot 
showed that RA treatment of wild type mice caused a 
downward shift in C2 from the untreated control group 
(Figure 3A). In contrast, an upward shift was found due to 
hepatic RXRa deficiency in comparison to the untreated 
wild type mice. Thus, RA treatment and RXRa deficiency 
had opposite effects. In addition, no significant change 
was noted when RA was used to treat hepatic RXRa KO 
mice. These findings unequivocally prove that the effects 



of RA on regulating those lipid homeostasis genes were 
RXRa dependent. Score plot (Figure 3A) indicates compo- 
nent 2 made a contribution to distinguish groups of con- 
trol, RA-treated, and RXRa-deficient mice. Thus, 114 out 
of 579 genes with high loading values (>0.5 or < -0.5) in 
component 2 were selected for further analysis (Figure 3B). 
Among them, 55 genes were induced by RA and had de- 
creased expression levels due to RXRa deficiency. The 
other 59 genes, whose expression levels were suppressed 
by RA, had increased expression levels due to a lack of 
RXRa. Thus, the expression levels of those 114 lipid- 
related genes are ligand (RA)-responsive and receptor 
(RXRa) -dependent. Based on the known function of those 
genes described in KEGG and PubMed, the role of those 
114 genes was assigned and summarized in Table 2. Re- 
markably, RA regulated many genes involved in certain 
pathways. For example, RA decreased the expression of 10 
genes in the cholesterol biosynthesis pathway, but did not 
increase the expression of any other genes in the same 
pathway. Thus, it is very likely that RA inhibited the bio- 
synthesis of cholesterol in an RXRa-dependent manner. 
RA also induced the expression of 13 genes in the RA 
elimination process and yet did not reduce the expression 
of any gene in the same process. Thus, RA can self- 
regulate its own level. Furthermore, RA also induced the 
expression of 9 genes in the biosynthesis of unsaturated 
fatty acids responsible for anti-inflammation. Since there 
was no inhibition of gene expression in the same pathway, 
it is very likely that RA up-regulates the synthesis of unsat- 
urated fatty acids and has an anti-inflammatory role. 
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Figure 3 PCA of the expression level of 579 lipid genes in wild type and hepatic RXRa-KO mice treated with and without RA. Wild type 
and hepatic RXRa KO mice were treated with and without RA for 7 days (1 50 mg/kg diet, n = 3). The expression of 579 genes involved lipid 
homeostasis was studied. (A) Score plot of PCA showing the difference among the groups. Spots within an ellipse belong to the same group. 
Arrows represent the direction separating groups from the wild type (WT) mice (open square □). There is no significance in the C1 component 
observed across all groups. (B) A histogram that shows the loading value of genes on C2. Genes with high loading value (> 0.5) on C2 had 
higher mRNA levels in RXRa KO than wild type livers. In contrast, genes with low loading value (< -0.5) on C2 had increased mRNA levels due to 
RA treatment. 
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Table 2 Biological functions responding to RA treatment 
and RXRa knockout in wild type mice 

Gene number 
(with RXRa bindings) 

Biological functions RA induced RA repressed 

& RXRa KO & RXRa KO 
repressed induced 



ChlP-Seq data indicated that most of the genes (87 out of 
1 14) had RXRa binding implying direct gene regulation. 

Taken together, ligand (RA) treatment and hepatic 
RXRa deficiency resulted in opposite effects. Figure 4 
summarizes the effect of RA and hepatic RXRa deficiency 
on lipid homeostasis. RXRa deficiency tends to favor satu- 
rated fatty acids, triglyceride, cholesterol, and bile acids 
synthesis. In contrast, RA treatment leads to unsaturated 
fatty acids and phospholipid synthesis and lipolysis as well 
as triglyceride breakdown. 

Binding of RA/RXRa responsive genes by other nuclear 
receptors 

Additional analysis was done to understand which other 
nuclear receptors may be involved in regulating the ex- 
pression of these 114 RA/RXRa target genes, which have 
a role in lipid homeostasis. The binding data generated in 
the current study (RXRa and RARa) were compared with 
the binding data of PXR, LXR, FXR, and PPARa. Figure 5 
shows overlapping genes with RXRa-heterodimers, as 
assessed by overlapping binding of RXRa and other nu- 
clear receptors. The data were organized by the number of 
different nuclear receptors binding the genes. For example, 
motifs located in the Abcal, Abhd5, Acsl, and Aldh3a2 
genes could be bound by RXRa and all five nuclear recep- 
tors. Peaks located in the Apoa4, Cyp51, Cyp7bl, and Elovll 
could be bound by RXRa and any 4 out of the 5 studied nu- 
clear receptors (Figure 5). Some of the commonly regulated 
genes have nuclear receptor binding site at the same loca- 
tion. The data indicated extensive crosstalk among nuclear 
receptors in regulating the expression of those genes. 

Quantification of serum cholesterol, triglyceride, and bile 
acid levels 

ChlP-Seq and RNA expression profiling indicate the role 
of RA in controlling lipid homeostasis in the liver. Serum 
cholesterol, triglyceride, and bile acid levels were quanti- 
fied to test the genetic findings. The data showed that RA 
reduced serum cholesterol, triglyceride as well as bile acid 
levels in wild type mice (Figure 6). However, such effects 
were not found in hepatic RXRa KO mice. In addition, 
serum cholesterol and triglyceride, but not bile acid, levels 
were elevated due to RXRa deficiency. These biochemical 
findings confirm the role of RA in regulating lipid homeo- 
stasis in the liver. 

Discussion 

This study establishes the role of nuclear receptors and RA 
in regulating lipid homeostasis in the liver. In addition, the 
mechanisms by which nuclear receptors and RA regulate 
lipid homeostasis were illustrated at the gene, pathway, and 
systemic levels. Furthermore, relationships between RXRa 
and PXR, LXR, FXR, as well as PPARa in regulating lipid 
homeostasis were analyzed. These nuclear receptors depend 



Lipid droplet growth 1 (0) 2 (2) 

Transportation of bile aicds for bile 1 (0) 0 
excretion 

Dehydrogenation of saturated fatty acids 3 (2) 0 
to unsaturated fatty acids 

Biosynthesis of glycerol phosphalipids 4 (3) 1 (1) 

Tryglyceride degration 4 (3) 0 

Biosynthesis of unsaturated fatty acids 9 (7) 0 
responsible for aniti-inflammation 

Elimination of retinoic acids 13 (11) 0 

Biosynthesis of bile acids 4 (4) 4 (4) 

Biosynthesis of retinoic acids 1 (1) 2 (2) 

Biosynthesis of saturated fatty acids 5 (5) 10 (10) 

Degradation of glycerol phosphalipids 1 (1) 4 (4) 

Degradation of saturated fatty acids 10 (10) 5 (5) 

Elimination of steroid hormones 3 (3) 5 (5) 

Fat digestion and absorption 3 (3) 7 (7) 

Lipid droplet breakdown 1 (1) 1 (1) 
(fat mobilization) 

Recycle of bile acids via hepatic-intestine 1 (1) 3 (3) 

S1P degradation 1 (1) 1 (1) 

Transportation of bile aicds for kidney 2 (2) 3 (3) 
excretion 

DHS1P degradation 1 (1) 0 

Elimination of unsaturated fatty acids 1 (1) 0 
(PGE2) responsible for lipolysis inhibition 

phosphatidylcholine to 1 (1) 0 
phosphatidylethanolamine 

Sphingolipid biosynthesis 1 (1) 0 

SPH (SM) degradation 0 1 (0) 

Biosynthesis of cholesterol 0 10 (10) 

(D 
(D 
(D 



Biosynthesis of steroid hormone 0 

Biosynthesis of tryglycerides 0 

Biosynthesis of unsaturated fatty acids 0 
responsible for pro-inflammation 

Breakdown of phosphalipid to form 0 
unsaturated fatty acids 

Elimination of cholesterol (from 0 
cyculation back to liver for catabolism) 

Elimination of cholesterol via steoid 0 
hormone pathway 

Phosphatidylethanolamine to 0 
phosphatidylcholine 



He ef al. BMC Genomics 2013, 14:575 
http://www.biomedcentral.com/1471-2164/14/575 



Page 7 of 1 1 



Oxidation 



RXRa KO induced 
CZ^ RA induced 

Elimination 



Lipolys 




Anti-inflammation 

Pro-inflammation 



Breakdown 



Elimination 



Breakdown 



Bile secretion 



Kidney secretion 
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on RXRa to execute their functions and more than 85% of 
their bound genes overlap with those bound by RXRa. The 
nuclear receptor binding data were strengthened by profil- 
ing the expression patterns of 576 lipid genes, which 
showed that RA treatment and RXRa-deficiency had an op- 
posite effect in regulating lipid homeostasis. Nuclear recep- 
tor binding data also illustrated extensive cross talk among 
the studied nuclear receptors. Thus, our novel in vivo data 
provided extensive evidence showing the role of RA in dic- 
tating lipid homeostasis in the liver. 

It is intriguing that more than 85% of the PXR, LXR, 
FXR, and PPARa binding sites overlapped with those of 
RXRa. In contrast, only 43% of RARa binding sites over- 
lapped with those of RXRa. It has been shown that RARa 
can form homodimers [22]. It is also possible that RARa 
may dimerize with RXR|3 and RXRy to exert its function. 
Although the percentage of overlapping bindings between 
RARa and RXRa is not as high as others, the number of 
genes that could be bound by RXRa/RARa (4554) is the 
highest, followed by RXRa/PPARa (3468), RXRa/FXR 
(2019), RXRa/LXR (988), and RXRa/PXR (666), implying 
the relatively extensive role of these nuclear receptors in 
regulating hepatic gene expression. It is important to fur- 
ther study the role of RARa and other RARs in the liver. By 
forming partners with other nuclear receptors, RXRa is a 
master regulator. Our data showed that more than 8000 
hepatic genes were bound by RXRa, and 72% of them over- 
lapped with the genes bound by RARa, PXR, LXR, FXR, or 
PPARa. The remaining 28% of RXRa binding sites might 
be bound by RXRa homodimer or the heterodimer of 



RXRa and VDR or CAR. Thus, the five nuclear receptors 
(RARa, PXR, LXR, FXR, or PPARa) analyzed in the current 
study account for almost three quarters of RXRa binding 
genes in the liver. Furthermore, nearly 50% of RXRa bind- 
ings overlapped with the bindings of PXR, LXR, FXR, and 
PPARa (Figure IB). Hence, lipid regulation should be one 
of the major functions of RXRa. 

Clustering and PCA showed that the genome-wide 
binding pattern of RARa is not similar to that of PPARa, 
LXR, PXR, and FXR. Biological function annotation also 
showed that RARa has some unique features including 
protein processing, protein localization, and RNA process- 
ing. However, the five studied nuclear receptors also dem- 
onstrated functional redundancy. For example, there are 
four pathways, including oxidation reduction, carboxylic 
acid catabolic process, organic acid catabolic process, and 
cofactor metabolic process, that can be regulated by more 
than four nuclear receptors. This finding suggests the im- 
portance of these four pathways in the liver, and the role 
of RARa in them. 

Although PPARa, LXR, FXR, and PXR have extensive 
roles in regulating lipids, they also have specific roles in 
regulating different types of lipids. RXRa/PPARa prefers to 
bind to genes that participate in neutral lipids, glycerol 
ether, and organic ether as well as fatty acid metabolism 
processes. All of which are either fatty acid-derived prod- 
ucts or precursors for the biosynthesis of fatty acids. An- 
other pathway bound by RXRa/PPARa is the acylglycerol 
metabolic process, which is involved in triglyceride homeo- 
stasis. RXRa/LXR tends to regulate genes involved in sterol 
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RXRa-RARa 



RXRa-PXR 



RXRa-LXR 



RXRa-FXR 



RXRa-PPARa 



Bound by RXRa and other 5 
NRs (30 genes): Abcal, Abhd5, 
Acsll, Aldh3a2, Aldh9a1, 
Angptl4, Aqp4, Cbr1, Cds2, 
Cptla, Cyp26a1, Cyp8b1. Dgat2, 
Ehhadh, Elovl2, Fabpl, Hsd11b1, 
Hsd17b4, Ldlr, Nr0b2, Plin5, 
Ppap2b, Ppara, Retsat, Scp2, 
Slc22a1,Slc27a2, Slco1a1, 
Tm7sf2, Ugt2b34 



Bound by RXRa and other 4 
NRs (49 genes): Apoa4, Cyp51, 
Cyp7b1, ElovM, Fabp5, Fdftl, 
Glyctk, Got2, Gpd2, Gpx1, 
Lpcat3, Neu1, Pnpla7, Rdh11, 
Rxra, Sqle, Sult1e1, Tecr, Ubc, 
Ugt2b1, Ugt2b5, Cyp2b10, 
Cyp2c29, Cyp2c37, Cyp2c50, 
Atp1b1, Scd2, Abcc3, Pla2g6, 
Ppap2c. Abcb4, Acaala, Acaalb, 
Acadm, Acot8, Adcy9, Aldh7a1. 
Aqp8, Aqp9, Chptl, Elovl6, 
Etnkl, Fabp2, Hsd17b12, 
Hsd17b7, Mgll, Pemt, Prkaca, 
Sc5d 



Bound by RXRa and other 3 
NRs (12 genes): Cyp2c70, 
Slco1b2, Dhrs4, Hadhb, Nr1h3, 
Rdh16, Sgppl, Cyp2j5, Cyp7a1, 
Scd1, Apoa5, Lipa 



Bound by RXRa and other 2 
NRs (3 genes): Cyp4f14, Acadl, 
Apoc3 



Uniquely bound by only one NR 
(2 genes): Amacr, Slc22a7 




Figure 5 Heat maps showing the binding profiles of RXRa/nuclear receptor on 96 RA-responsive and RXRa-dependent lipid genes. RA 

responsive and RXRa-dependent genes (1 14) were analyzed by nuclear receptor binding. Among those 1 14 genes, 96 genes shown in black lines 
can be commonly regulated by RXRa and its partner. Some of the commonly regulated genes have nuclear receptor binding site at the 
same location. 
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metabolism, which is consistent with its known role [23]. 
RXRoc/FXR not only binds to the genes participating in 
steroid metabolism process, but also those involved in lipid 
transport and carbohydrate metabolism processes. RXRa/ 
PXR binds to the genes involved in regulating the pyruvate 
metabolic process at the DNA binding level. Pyruvate is a 
key intersection for fatty acid, carbohydrate, and protein 
metabolisms. In addition, RXRoc/PXR also regulates re- 
sponse to acute phase, inflammatory, and wounding, im- 
plying that PXR can be an excellent target for metabolism 
and inflammation-related health issues. Lastly, FXR binds 
to the genes involved in monosaccharide metabolism, 
which shows the intimate relationship between bile acid 
and glucose homeostasis [24]. 

All of the 114 genes that showed differential effects 
of RA treatment and RXRa deficiency are bounded 
by RXRa and RARa, PXR, LXR, FXR, and PPARa 
heterodimers. These findings indicate that those stud- 
ied nuclear receptors retain RA response in vivo and 
the effect of RA is dependent upon those nuclear re- 
ceptors. RA has a broad spectrum of effects including 
biosynthesis of retinoids, phospholipids, and unsatur- 
ated fatty acids. It also has a role in eliminating reti- 
noids, oxidizing saturated fatty acids, and breaking 
down triglycerides. It seems that RA has extensive 
beneficial effects in maintaining the health of the liver. 
Specifically, RA induced the expression of Cyp2c37/ 
38/50/54/70 and Cyp2j5. These genes encode enzymes 
involved in the generation of epoxyeicosatrienoic acids 
[25], which have anti-inflammatory effects [26]. In 
contrast, RXRa deficiency induces the gene expression 
of Cyp4f that is responsible for the generation of 20- 
hydroxyeicosatetraenoic acid, a pro-inflammation mol- 
ecule [25]. In addition, RA increases mRNA levels of 
cbrl(carbonyl reductase 1), which is responsible for 
transforming prostaglandin E2 to prostaglandin F2a. 
Prostaglandin E2 and F2a have different effects in 
regulating lipid breakdown. Prostaglandin E2 is a lip- 
olysis inhibitor [27]; whereas, prostaglandin F2a has 
not been shown to have the same effect. Thus, the in- 
duction of cbrl gene expression could be a mechanism 
by which RA induces lipolysis. RA also induces expres- 
sion levels of gene encoding proteins for phospholipid 
biosynthesis, but RXRa deficiency increases the expres- 
sion of the genes that have a role in the degradation of 
phospholipids. This finding suggests the potential role 
of RA in maintaining the normal structure of the cell 
membrane. Formation of the monolayer of lipoprotein 
or lipid droplet is one of the major ways that phospho- 
lipids regulate lipid metabolism [28]. Phosphatidyletha- 
nolamine (PE) and phosphatidylcholine (PC) are two 
important phospholipids that show different effects on 
lipid metabolism in humans and rodents. Lower PC/PE 
ratio induces steatosis or even steatohepatitis in humans 



[29], however, PE has a greater effect than PC in redu- 
cing the cholesterol level in rodents [30]. Our data 
showed that RA induced the gene expression of ptdssl 
(phosphatidylserine synthase 1), which converts PC to 
PE. Consistently, the ptdssl gene expression is reduced 
due to hepatic RXRa deficiency. These gene expression 
levels and DNA binding data not only showed the under- 
lying mechanism for RA in regulating liver gene expres- 
sion, but also suggested the biochemical outcome. 

Conclusions 

Taken together, the differentiation and morphogenetic 
effect of RA is well known. However, the current study 
provides a comprehensive analysis of the role of RA in 
lipid homeostasis. All-trans RA is the most abundant 
retinoid that can be easily detected in the liver. Thus, 
RA is likely to act as a regulator to control hepatic lipid 
metabolism. Since the effect of RA is broad, it is import- 
ant to develop specific retinoids in order to target spe- 
cific pathways. Such efforts may allow us to identify 
compounds that can be used to treat or prevent meta- 
bolic syndromes and other lipid-related health issues. 

Methods 

Materials and animal models 

Male wild type mice (12 weeks old) and hepatocyte 
RXRa-deficient mice (KO) [10,11], which have the same 
genetic background of C57BL/6, were used. The RXRa 
KO mice were produced and characterized previously 
[10,11]. The LoxP sites were inserted into introns flanking 
the fourth exon of the RXRa gene covering the DNA 
binding domain, which is deleted after crossing the floxed 
RXRa allele against a transgenic line in which ere recom- 
binase is expressed under the control of the albumin pro- 
moter. The mutant mice express a truncated protein that 
has the intact ligand binding domain, but lacks the DNA 
binding domain. Animal protocols and procedures were 
approved by the Institutional Animal Care and Use Com- 
mittee (IACUC) at the University of Kansas Medical Cen- 
ter and the University of California, Davis. 

All ChlP-grade antibodies except anti-RNA Pol II 
(Millipore, MA) were purchased from Santa Cruz Biotech- 
nology, Inc. (Santa Cruz, CA). DNA purification kit was 
purchased from Qiagen Co. (Valencia, CA). All other ChlP- 
related reagents were obtained from Invitrogen Co. (Carls- 
bad, CA). 

Animal treatment 

Retinoic acid (Sigma-Aldrich, MO) was given to wild 
type and hepatic RXRa-deficient mice at a dosage of 
150 mg/kg diet for 7 days. As controls, mice were fed 
with normal diet. 
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Chromatin immunoprecipitation (ChIP) 

ChIP was performed according to our previously pub- 
lished study [31]. After fixation, the mouse livers were 
subjected to lysis with cell and nuclear lysis buffer. Sonic- 
ation was used to fragment the chromatin, followed by 
precipitation with specified antibodies. The target DNA 
fragments were obtained by reverse crosslinking and puri- 
fication. Antibodies against IgG and RNA Pol II were used 
as negative and positive controls, respectively. 

DNA library preparation and sequencing 

By using the End-It DNA End Repair Kit (Illumina, 
Madison, WI), DNA fragments prepared from ChIP 
were ligated with specified adaptors and amplified, then 
size-selected (175-225 bp) on an agarose gel followed by 
sequencing (High-Seq 2000, Illumina, Madison, WI). 

Alignment, call peak, and annotation of ChlP-seq data 

The target sequences were aligned to the mouse genome 
(http:/ /hgdownload.cse.ucsc.edu/goldenPath/ mmlO/bigZips/) 
by Bowtie 0.12.7 [32] followed by peak-calling using MACS 
(version 1.4.1) [33]. The peaks were annotated using the 
database (NCBI37/mm9) by Peak Analyzer [34]. The back- 
ground cut off standard was set to be 20 fold of the input sig- 
nals [18]. The cut off distance from the transcription start 
site (TSS) was set to be 10 kb. Co-localization is defined as 
having at least 25% overlap in their peak widths. 

Microarray 

Affymetrix 430 A_2 Chip (Santa Clara, CA) was used to 
determine the genome-wide mRNA expression levels. 
Microarray data were annotated using Affymetrix Expres- 
sion Console (MAS5). The probe signal with p values less 
than 0.05 were used for further analysis. 

ChlP-seq data analysis 

All data were treated with the same cut off criteria. The 
generated RXRoc binding data were compared with the data 
for RARoc, PXR, LXR, FXR, and PPARa. The principle 
component analysis (PCA) and cluster analysis package in 
SPSS program was used to analyze the global binding data. 
For both PCA and cluster analysis, called peaks were 
assigned the value of 1. Not called peaks were assigned the 
value of 0. Genes with overlapping binding sites of RXRa 
and each of RARa, PXR, LXR, FXR, and PPARa at the 
same location were functionally analyzed by the DAVID 
(http:/ /david.abcc.ncifcrf.gov/) [35] . 

Lipid homeostasis analysis based on mRNA expression 

Genes (579) involved in regulating lipid homeostasis were 
extracted from the KEGG database (Kyoto Encyclopedia of 
Genes and Genomes, http://www.genome.jp/kegg/). The 
expression of those 579 genes were determined in wild 
type and liver RXRa KO mice treated with and without 



RA (n = 3) for 7 days. After multiple comparisons, only 30 
and 36 out of all 579 lipid homeostasis genes showed sig- 
nificant change at the corrected p-value of 0.05 after RA 
treatment and RXRa knockout, respectively. Therefore, 
IBM SPSS PCA package was used to differentiate groups 
based on the global expression pattern of all 579 lipid 
homeostasis genes. 

Serum lipid assays 

Triglyceride, cholesterol, and bile acids in the serum were 
assayed using a commercially available kit (Pointe Detroit, 
Michigan) that was modified to a 96-well format. Spectro- 
photometric analysis was conducted with a Bio-Tek mi- 
cro titer plate reader (Bio-Tek, VT). 

Availability of supporting data 

The microarray and ChlP-Seq data supporting the re- 
sults of this article are available with accession numbers 
of GSE50028 and GSE46762, respectively, in the GEO 
repository (http://www.ncbi.nlm.nih.gov/geo/). 
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